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ABSTRACT 

The recently discovered Virgo stellar over-density, which expands over ~1000deg 2 perpendicularly 
to the Galactic disk plane (7< Z < 15 kpc, R^ 7 kpc), is the largest clump of tidal debris ever 
detected in the outer halo and is likely related with the accretion of a nearby dwarf galaxy by the 
Milky Way. We carry out N-body simulations of the Sagittarius stream to show that this giant stellar 
over-density is a confirmation of theoretical model predictions for the leading tail of the Sagittarius 
stream to cross the Milky Way plane in the Solar neighborhood. Radial velocity measurements are 
needed to confirm this association and to further constrain the shape of the Milky Way dark matter 
halo through a new generation of theoretical models. If the identification of Virgo over-density and 
the Sagittarius leading arm is correct, we predict highly negative radial velocities for the stars of Virgo 
over-density. The detection of this new portion of the Sagittarius tidal stream would represent an 
excellent target for the on-going and future kinematic surveys and for dark matter direct detection 
experiments in the proximity of the Sun. 

Subject headings: galaxies: individual (Sagittarius) — galaxies:interactions — Galaxy: structure - 
Galaxy: halo 



1. INTRODUCTION 

In the last decade, various large scale CCD surveys 
(Sloan Digital Sky Survey, 2MASS 7 , QUEST) have for 
the first time found strong evidence for the presence of 
a significant amount of substructure in the halo of the 
Milky Way in form of long tidal streams or stellar clumps 
(Newberg et al. 2002; Majewski et al. 2003). These have 
been interpreted as the remnants of latest mergers occur- 
ring as a part of the hierarchical build-up of our Galaxy. 
In addition, N-body simulations of these merger events 
(Helmi 2004; Martfnez-Delgado et al. 2004a; Law, John- 
ston & Majewski 2005; Penarrubia et al. 2006) have 
been extensively used to, firstly, associate the tidal de- 
bris with their progenitor galaxies, secondly, determine 
the dynamical history of the progenitor satellites and, 
thirdly, constrain the distribution of dark matter in the 
Milky Way. The ultimate goal of these studies is to infer 
the formation process of the Milky Way from the prop- 
erties of the present fossil records. 

Much of the observational and theoretical effort has 
focused specially on the two largest tidal streams dis- 
covered so far: the tidal stream of the Sagittarius (Sgr) 
dwarf galaxy, which wraps around the Galaxy in a highly 
inclined orbit with respect to the Milky Way disk (Ibata 
et al. 2001; Martmez-Delgado et al. 2001; Majewski 
et al. 2003) and the Monoceros tidal stream (Yanny et 
al. 2003; Ibata et al. 2003), a low-latitude tidal stream 
whose progenitor, possibly the controversial Canis Major 
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dwarf galaxy (Martin et al. 2004; Bellazzini et al. 2004; 
Martfnez-Delgado et al. 2005), moves on a prograde, 
nearly circular orbit (Penarrubia et al. 2005), according 
to the recent proper motions measured by Dinescu et al. 
(2005). Additionally, Rocha-Pinto et al.(2004) reported 
the discovery of a new possible tidal stream in the Trian- 
gulum/Andromeda constellations, although subsequent 
theoretical simulations identified this tidal debris as a 
part of a more distant, metal poor wrap of the Mono- 
ceros tidal stream progenitor (Penarrubia et al. 2005; 
Martfnez-Delgado et al. 2007, in preparation). 

In addition to these two large tidal streams, the 
analysis of large samples of stellar tracers of halo sub- 
structure (RR Lyrae, F-type turnoff stars) have revealed 
the presence of numerous stellar over-densities possibly 
associated with unknown merger events in the outer 
halo. In the Virgo constellation, several groups have 
reported numerous detections of stellar over-densities 
using different tracers: 

(i) The most conspicuous one is a diffuse concentration 
of metal poor, old stars situated at ~ 20 kpc from the 
Sun in the direction of the Virgo constellation, first 
detected by means of the identification of its associated 
main-sequence turnoff (MSTO) in a color-magnitude 
diagram (CMD) of a SDSS slice (named S297+63-20.0; 
Newberg et al. 2002). Additionally, Vivas & Zinn (2003) 
reported a clump of 21 RR Lyrae stars located at ~ 
20 kpc from the sun in the same region of the sky (see 
also Ivezic et al. 2005 and Vivas & Zinn 2006). This 
over-density (named the 12 h .4 clump) extends over an 
area of the sky of 15° x 1° (corresponding to to a spatial 
size of 5.3 kpcx0.3 kpc at D ~ 20 kpc) and display a 
line-of-sight depth of only ~ 1.5kpc (Vivas 2002), which 
is considerably smaller than that of a similar RR Lyrae 
over-density detected at the apocenter of the Sgr stream 
(~ 4 kpc; Vivas 2002). Owing to its possible tidal origin, 
they denoted it as the Virgo stellar stream (VSS). 
Martfnez-Delgado et al. (2004a) suggested that the 
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VSS might be related to the complex debris structure 
of the Sgr stream predicted by the theoretical models in 
this region of the sky. However, spectroscopic follow-up 
of a sample of these RR Lyrae stars by Duffau et al. 
(2006) showed a fairly narrow velocity distribution that 
pointed to a progenitor galaxy less massive than the 
Sgr dwarf and, therefore, to a hitherto unknown tidally 
disrupted dwarf galaxy. Using SDSS photometric data, 
these authors estimated that the VSS covers roughly an 
area of 106 deg 2 in the sky. 

(ii) More recently, Juric et al. (2006) detected 
a remarkable excess of main-sequence stars towards 
the direction of Virgo constellation, named the Virgo 
over-density (VOD). These authors construct a three- 
dimensional density distribution of Milky Way stars 
based on the position of ~ 48 million stars obtained by 
the SDSS in order to identify stellar over-densities. The 
VOD density map reveals a giant, diffuse lump of stars 
without any evidence of a distinct core that extends over 
more than ~ 1000 deg 2 (see Fig. 1) towards the direc- 
tion (l,b) — (300°, 65°). This stellar over-density shows 
an remarkably extended vertical structure perpendicular 
to the Galactic plane (AZ ~ 8 kpc, with heliocentric dis- 
tances between 5-17 kpc) that possibly crosses the Galac- 
tic plane, extending into the southern Galactic hemi- 
sphere (see also Fuchs, Phleps & Meisenheimer 2006). 
The density of VOD main-sequence stars peaks at ~ 16 
kpc from the Sun (or RR Lyrae magnitude of ~ 16.7), in 
good agreement with the RR Lyrae distance determina- 
tion. These authors suggest that the giant stellar over- 
density might be related with a hitherto unknown merger 
event of a nearby, low metallicity dwarf galaxy .Alterna- 
tively, Xu, Deng & Hu (2006) suggest that this stellar 
lump could be related to a non-axisymmetric Galactic 
stellar spheroid (see Sec. 13. lj) . 

Although VOD and VSS are coincident in position and 
display a huge sky-projected size, the different vertical 
structure suggests that they are independent systems. In 
absence of radial velocity measurements for the VOD, it 
is not possible to address whether they share a common 
origin or are completely unrelated. For this reason, we 
will refer to them with different names in this paper, 
following the nomenclature given to each over-density in 
the original papers. 

The main aim of this paper is to investigate the nature 
of the VOD and to analyze whether this stellar lump 
could be related to the Sgr stream through a comparison 
between up-to-date theoretical models and the available 
observational data. The structure, kinematics and pos- 
sible origin of the VSS will be analyzed in a companion 
paper. 

In Sec. 2 the theoretical models of Sgr tidal stream 
are described. In Sec. 3 we discuss the possible origins 
of the VOD. In particular, we explore the presence of 
tidal debris from the Sgr tidal stream in Virgo by means 
of update theoretical models and compare its expected 
structural parameter with those derived for Virgo. Fi- 
nally, conclusions are given in Sec. 4. 

2. THEORETICAL MODELS 

The giant stellar over-density detected in the Virgo 
constellation is located in a region of the Milky Way 
where several theoretical models predict the presence of 



Sgr galaxy debris. In order to explore a possible associa- 
tion of the VSS and VOD with the Sgr tidal stream, we 
shall compare their derived spatial and kinematic prop- 
erties with those obtained from up-to-date theoretical 
models. 

Our Galaxy model consists of a Miyamoto-Nagai 
(1975) disk, a Hernquist (1990) bulge and a Navarro, 
Frenk & White (1996) dark matter halo (hereafter, NFW 
halo) . The gravitational potential of each of those com- 
ponents in cylindrical coordinates is: 
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There, r v j r ,r s are the virial and scale radii, respectively, 
M h = M h (r vir ) and r 2 = R 2 + z 2 . The potential of 
an axi-symmetric NFW halo was calculated from Chan- 
drasekhar (1960) using elliptic coordinates 
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where qh is the axis-ratio of iso-density surfaces. In this 
contribution, we shall denote oblate and prolate halos as 
those with qh < 1 and qh > 1 , respectively. For spherical 
halos Eq. ([3]) reduces to 
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Following Johnston et al. (1999) we fix the disk and 
bulge parameters as M d = 1.0 x 10 n M Q , M b = 3.4 x 
10 10 M o , a = 6.5kpc, b = 0.26 kpc and c = 0.7kpc. The 
Milky Way halo parameters at z = were taken from 
Klypin, Zhao & Somerville (2002) being M h = 1.0 x 
10 12 Mq, r v ; r = 258kpc, r s = 21.5kpc, which leads to a 
concentration at the present epoch of c = r v - u /r s = 12. 
According to the results of Peharrubia et al. (2006), the 
evolution of the host Galaxy potential is not reflected in 
the present properties of tidal streams, so that we use a 
static Milky Way potential for simplicity. 

The actual disruption of satellites is modelled by 'live' 
(i.e. self-consistent, self-gravitating) N-body realiza- 
tions of a King (1960) model, with a dimensionless 
central potential Wq — 4, or concentration parameter 
c = log 10 (r t /rfe) ~ 0.84, where and r t are the King and 
tidal radii, respectively. Our satellite models have N = 
10 5 particles. The initial and final (i.e present) masses 
are M s (t ) = 10 9 M© and M s (t f ) = 5 x 10 8 M Q , respec- 
tively. The initial King and tidal radius are fit (to) = 0.58 
kpc and r t (to) = 4.01 kpc. 

The equation of motion for each satellite particle is: 



d 2 ^ 
dt 2 



(6) 
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where $ s is the self-gravitational potential of the satellite 
galaxy and i = 1, N. We use superbox (Fellhauer et 
al. 2000) to calculate <!> s at each time-step and solve Eq.[5] 
through a leap-frog scheme with a constant time-step of 
At = 0.65 Myr, which is about l/100th the satellite's 
dynamical time. 

Note that the Galactic potential is static and, therefore, 
Eq.[6]does not implement the response of the Milky Way 
to the presence of Sgr. We also neglect the effects of 
dynamical friction on the Sgr's orbit (as suggested by 
Law et al. 2005). 

We have reproduced the N-body models presented in 
Law et al. (2005). These models reproduce the following 
observational constraints 8 : 

1. The heliocentric position and radial velocity of the 
Sgr dwarf are {D,l,b) = (24kpc, 5.6°, -14.2°) and 
v r = 171 km/s, respectively. 

2. The orbital plane of Sgr has an inclination with 
respect to the Milky Way disk of i ~ 76°. 

3. The averaged heliocentric distance for Sgr leading 
debris is D ~ 50 kpc. 

Items 1. and 2. above show that we have observational 
measurements for five of the six coordinates that deter- 
mine the Sgr orbit. Law et al. (2005) surveyed the un- 
known coordinate (namely, the tangential velocity com- 
ponent, Dtan) so that the resulting N-body models repro- 
duced all the observational constraints listed. 

Among other free parameters, Law et al. (2005) also 
explored which halo axis-ratio (qu) would produce the 
best-fitting model to the available observational data. 
They found that, whereas the geometry and kinematics 
of the trailing arm are scarcely sensitive to the adopted 
qh , radial velocity measurements of the leading arm were 
best matched by prolate (qh > 1) halo models, in agree- 
ment with Hclmi (2004). However, Johnston, Law and 
Majewski (2005) showed that Sgr models in prolate ha- 
los cannot reproduce the observed precession rate in 
the youngest pieces of the Sgr stream. Furthermore, 
they were able to constrain the halo axis-ratio to be 

qh = 0.83 0.92, excluding models with qh > 1 at 

a 3<T level. In view of these results, we have performed 
N-body simulations where the halo axis-ratio was fixed 
either to qh = 0.8 or to qh = 1-4, excluding for the sake of 
brevity spherical halos, which neither reproduce the ra- 
dial velocity variation along the stream nor its precession 
rate. 

In a coordinate system where the Sun is located at 
(X, Y, Z) = (8, 0, 0) kpc, with a velocity of [U, V, W) = 
(-10, -220, 7) km/s (Binney & Merrifield 1998), the Sgr 
velocities that best fit the kinematic and spatial distri- 
butions of stream debris are (237,-35,220) km/s for 
q h = 0.8 (oblate halo) and (244,-39,249) km/s for 
qh = 1.4 (prolate halo). 

With those inputs, we have evolved the N-body satel- 
lite model from 4 Gyr in the past to the present. 

8 It is also important to remark that we have not input any 
observational data of the VOD in the fitting of these simulations, 
that are only based in the fitting of an independent set of data of 
the Sgr tidal stream. 



3. THE ORIGIN OF THE VIRGO OVER-DENSITY 

In this Section we analyse the likelihood of various sce- 
narios that may shed light on the nature of the VOD. 

3.1. The signature of a tri-axial Galactic stellar halo 

An alternative explanation for the presence of the un- 
expected MS feature in the CMD reported in several 
studies (and so-called the VOD: Newberg et al. 2002; 
Duffau et al. 2006; Juric et al. 2006; see Sec.l) is to 
postulate the existence of a non-axisymmetry Galactic 
model component, such as a triaxial halo (Newberg & 
Yanny 2005; Xu et al. 2006). We have explored this 
possibility by inspecting a sample of deep CMDs (with 
limiting magnitude Ru m ~ 24.5) of Galactic fields taken 
with the Wide Field Camera (WFC) at the prime focus of 
Isaac Newton 2.5-meter telescope (La Palma, Spain) dur- 
ing different observing runs devoted to the search of tidal 
debris associated to known Galactic satellites or globular 
clusters (Martfnez-Delgado et al. 2004b). Table 1 lists 
the position of our fields with positive and negative de- 
tections of the MS feature in the CMD. These negative 
detections can be interpreted as the lack of tidal debris 
or that its surface brightness is extremely low to be de- 
tected in the small field of view of the WFC (35' x 35') 9 . 

Fig. Q] shows the position of these negative detection 
(open circles) over-plotted in the density map of VOD 
derived by Juric et al (2006) from SDSS observations. 
The most interesting result is a negative detection ~ 20° 
North of the center of our target field (see Table 1) that 
clearly shows a decrease of stream star counts in a small 
area of the sky. Although our sparse pencil beam survey 
cannot constrain the spatial extent of this stellar popu- 
lation in the sky, the distribution of negative detections 
strongly indicates that this MS population is not likely to 
be a manifestation of a non-axisymmetric Galactic model 
component, but is probably associated with the remnants 
of a dwarf galaxy under tidal disruption, in agreement 
with previous studies. Juric et al. (2006) have reached 
the same conclusion by analyzing the vertical (Z) distri- 
bution of VOD debris from their sample of SDSS G stars 
(0.10 < r-i < 0.15). 

3.2. A tidal debris of the Sgr dwarf? 
3.2.1. Sky projected density maps 

It is well understood both from observations and the- 
oretical modeling that the Sgr leading arm turns around 
and falls towards the Galactic plane in the proximity of 
the Solar neighborhood, yielding a vertical debris struc- 
ture that might cross the Galactic plane and possibly 
extend into the southern Galactic hemisphere (Majew- 
ski et al. 2003; Martinez-Delgado et al. 2004a; Law et 
al. 2005). This downward flow of the Sgr stars onto the 
Galactic disk may share similar characteristic to those 
derived for the VOD in the stripe covered by the SDSS. 
As a guide for the reader, we illustrate this putative sce- 
nario in Fig. [2] together with the identification of the 
different parts of the Sgr tidal stream obtained from our 
theoretical models (see Sec. [5] for a three-dimensional 
comparison) . 

9 Martinez-Delgado et al. 2004a estimated the limit SB detection 
of our CMD technique is S ~ 32 mag arcsec -2 
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To examine whether the VOD is a part of the Sgr tidal 
structure, we compare the predictions from the theoreti- 
cal models outlined in Sec.Oto the SDSS data presented 
in Juric et al. (2006). When comparing the theoretical 
stream models against the SDSS, it is important to keep 
in mind that the observational data may be biased by 
systematic over or underestimates of the distance scale. 
The stellar distances in SDSS data were derived using 
an assumed photometric paralax relation whose shape is 
well constrained but the zero-point may be offset by as 
much as 0.4mag (Juric et al. 2006). This would translate 
into 20observed center of VOD to anywhere from 8 to 12 
kpc. However, even the worst case scenario would be 
compatible with the conclusions that we are to present. 

Secondly, the vertical shape and density profile of VOD 
(Fig. 4, bottom left panel) depend strongly on the de- 
tails of the (still preliminary) Galactic stellar halo model. 
The density gradients along this vertical structure should 
therefore be interpreted with caution, especially at low Z 
(closer to the disk). There, the number density of Galac- 
tic stars dominates over the stars in VOD and even a 
few percent mismatch in the Galactic model may lead to 
an artificial increase (or decrease) of the residual stellar 
density ascribed to VOD. 

In Fig. [3] we plot two equal-area projections of stream 
particles from the models with oblate qh = 0.8 and pro- 
late qh — 1.4 halos (top and bottom, respectively). In 
these plots, the North Galactic pole is in the center, con- 
centric circles are lines of constant Galactic latitude b and 
the radial lines denote locations of constant Galactic lon- 
gitude I. The Galactic center is to the left. We restrict 
the theoretical sample to the distance range of Juric et 
al. (2006) data by including only particles at b > and 
5kpc < d < 15kpc. To estimate the two-dimensional 
projected density distribution, we convolve the particles 
with a 2D Gaussian a = 2.9 deg kernel. The resulting 
density is shown as varying shades of gray (black and 
white regions denoting the highest lowest density, respec- 
tively) and can be directly compared with the observed 
morphology of the VOD displayed in Fig [T] (that was 
made with the same scale). 

These maps show a conspicuous, diffuse lump of Sgr 
stream stars that extends over a few hundreds of square 
degrees, whose shape and extent resembles those of the 
stellar over-density displayed in Fig. [TJ The position of 
the over-density center (defined as the location of max- 
imum density) in the sky depends on the shape of the 
halo, and is consistent with simulations of Law et al. 
(2005), Fig. 13. The oblate model (q h =0.8) produces a 
better agreement with the SDSS data, with the over- 
density being centered at (l,b)=(277,63). This corre- 
sponds to a small offset, ~ 11 deg, with respect to the ob- 
served position of the VOD ((/, fc)=(300,65)) as reported 
by Juric et al. (2006). Taking into account the uncer- 
tainties in the Sgr stream modeling (see Sec. [2|) and in 
the SDSS data as noted above, the agreement in posi- 
tion is remarkable. This result also provides a natural 
explanation for the huge apparent extension of this tidal 
remnant, since we are observing through a stellar tidal 
stream that approaches the Sun in an almost radial di- 
rection in a localized area of the sky. 10 

10 An analogy can be found when a meteoroid stream crosses 
the Earth; see Nakamura et al. 2000. 



The prolate (q/ l =1.4) halo model (bottom panel) pro- 
duces a slightly worse estimate of the location of VOD 
center and its general shape. It does predicts a substan- 
tial amount of faint, diffuse debris in the 180 < I < 240 
region. This is to a degree observed in SDSS data - e.g. 
the density in 180 < I < 210 (see Fig. [1} is apparently 
greater than that in its symmetric counterpart with re- 
spect to the I — 0, 1 = 180 direction. 

Interestingly, Gondolo et al. (2005) estimated a lo- 
cal stellar density of Sgr leading tail of between 270-740 
stars/kpc 3 , corresponding to a range of 2.1 to 7.4 x 10~ 7 
stars/pc 3 or -15.4 < ln(S) < -14.1. This is within the 
measured limits of stellar density in the VOD (see Fig. 
23 in Juric et al. 2006). 

We note that the interpretation of the faint density 
excess is muddied by the existence of a different narrow 
stream in that region (Juric et al. 2006), apparently un- 
related to Sgr stream. An alternative scenario is that 
the narrow stream might be the true origin of the den- 
sity excess seen in SDSS data, instead of the Sgr debris. 
However, as it will become evident in the next Section, 
the three-dimensional structure of the VOD strongly dis- 
favors this hypothesis. 

3.2.2. Three-dimensional spatial structure 

We compare now the three-dimensional structure of 
the VOD to the structure obtained from the theoreti- 
cal models of the Sgr stream in that region of the sky. 
We examine various cross-sections (horizontal and per- 
pendicular to the Galactic plane) of the VOD and the 
equivalent cross-sections of the Sgr stream model. 

In Fig. 2] we plot the measured stellar density (left) and 
that obtained from the Sgr stream models for two differ- 
ent values of halo flattening, q h =0.8 (middle) and qh=lA 
(right). The density shown on the top left panel is that 
shown in Fig. 8 (middle-left panel) of Juric et al (2006). 
Before discussing its implications, we will very briefly de- 
scribe the method used here (for a more detailed discus- 
sion we refer the reader to Juric et. al 2006): The SDSS 
has observed ~ 4.2 million G stars (0.1 < r — i < 0.15) 
within the survey area. To convert the SDSS star counts 
to spatial number density, we estimated the distances to 
individual stars using the photometric parallax relation 
of Juric et al. (Eq. 1). From the estimated absolute and 
measured apparent magnitudes, and from astrometry, we 
calculate the three-dimensional (X, Y, Z) , positions of all 
stars in the SDSS stellar sample. We then bin them to 
dx = 500pc pixels in the coordinate system as defined 
above to obtain a three-dimensional star-count map. Fi- 
nally, by dividing the map by the volume that was ob- 
served within each voxel, we obtain a three dimensional 
map of number density within the Milky Way. 

The density shown in the left panel of Fig. [l]is a planar 
cross-section of the 3D map along a plane parallel with 
the Galactic plane, and at Z — lOkpc above it. The col- 
ors encode the stellar number density (stars pc~ 3 ) on a 
natural log scale, ranging from cxp(— 15.46) = 1.910~ 7 
(blue) to exp(-14.03) = 8.310" 7 stars pc~ 3 (red). The 
black regions are places where SDSS has not observed. 
The yellow dot marks the position of the Sun in X-Y 
plane (X — 8kpc, Y — Okpc) The circles visualize the 
presumed axial symmetry of the Galaxy. If the Galactic 
number density distribution was axially symmetric, the 
color contours in the observed regions would trace such 
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circles around the Galactic center, and the density dis- 
tribution would respect the symmetry with respect to re- 
flection along the X axis (p(X, Y, Z) = p(-X, Y, Z)). As 
seen in the left panel, the density distribution is neither 
circular nor symmetric, and exhibits a strong density en- 
hancement centered approximately at X — 6, Y — 4kpc, 
the VOD (Juric et al. 2006). 

We test under which conditions the VOD might be 
plausibly identified with the leading arm of the Sgr 
tidal stream by over-plotting the equivalent X-Y slices 
at Z = lOkpc of simulated Sgr debris over the contours 
of the SDSS observed area (the gray area in middle and 
left panels). The figures in the middle and left panels are 
constructed in a manner similar to the method used for 
SDSS data outlined above. However, we note that the 
particle counts have been averaged in the 9 kpc < Z < 
11 kpc slice. The difference in thickness, compared to 
dx = 500pc thickness used for SDSS density estimation is 
to account for the effect of radial smearing introduced by 
the intrinsic <jm t ~ 0.3mag dispersion of the photomet- 
ric parallax relation. This effect increases the effective 
radial SDSS voxel scale (for a more detailed discussion, 
see Juric et al. 2006). We estimate the density of Sgr 
debris by convolving the simulation particles in the ob- 
served volume with a a — 0.5 kpc Gaussian kernel, thus 
obtaining the final distribution shown in the panels. The 
color scheme is that used in the right panel, with equal- 
density contours added. We also plot, as small yellow 
crosses, the individual simulation particles that fall off of 
the lowest equal-density contour. 

As it was the case in Fig. [31 the oblate qh = 0.8 model 
agrees well with the data within the observational and 
theoretical uncertainties. It successfully reproduces the 
diffuse shape of the over-density and the location of over- 
density maximum. The location of the maximum for 
qh = 0.8 model locates within 1.5 to 2kpc of the ob- 
served R ~ 6.5 — 7kpc maximum of VOD, and the dif- 
fuse morphology of the qh = 0.8 model appears to show 
a considerably better match to the observed morphology 
than that of the qh — 1.4 model. Some debris from the 
qh = 1.4 model (right panels) can be also found at the 
location of the VOD. However, this model predicts the 
existence of other over-densities where none have been 
detected and the morphology at the VOD location is ir- 
reconcilable with that inferred from the SDSS data. 

We now examine the vertical distribution of matter in 
the VOD and compare it with the models of Sgr stream. 
The bottom panels of Fig. [4] compare the observed verti- 
cal stellar number density distribution in the Galaxy, cut 
along a narrow <f> — 30 deg plane that passes through the 
center of the Galaxy and the center of the VOD (the red 
dashed box on top middle and left panels of Fig. [2| . The 
horizontal coordinate is the distance from the Galactic 
center, R, along the 4> = 30° direction, while the vertical 
coordinate is the height above the Galactic plane. As 
in the case of the top left panel, the bottom left panel 
shows the over-density as measured from SDSS data (and 
is equivalent to the top right panel of Fig. 22 in Juric et 
al.2006), while the two panels to the right show the sim- 
ulated density of the Sgr tidal stream as obtained from 
our simulations (Sec. [2|) using an oblate (^=0.8) and a 
prolate (^=1.4) halo. The color scheme is the same as 
for the top panels, except that in the case of the bottom 
left panel we show the fractional residuals with respect 



to the best fit Galactic halo model of Juric et al. (2006). 
In the case of qh=0-8 halo, the Sgr leading arm is nearly 
perpendicular to the Galactic plane, with a slight down- 
ward density gradient, in agreement with the findings of 
Juric et al. (2006). Its location agrees to within AR ~ 1 
kpc of the observed position of Virgo over-density. In 
contrast, although the prolate models (g^ = 1.4) correctly 
reproduce the Galactocentric distance of the over-density 
(R r-j 6 — 7kpc), they fail in order to reproduce the ob- 
served morphology. These models show a strong den- 
sity maximum at (R, Z) ~ (7, 8)kpc (incompatible with 
the density map) and pieces of debris scattered out to 
R ~ 12.5kpc (not detected by SDSS). In this working 
scenario, the VOD morphology strongly disfavours pro- 
late halo models. 

The conclusions obtained from these cross-sections can 
be summarized in the different spatial perspectives dis- 
played in Figure^ 1 , which clearly shows the coincidence 
in position of the VOD and the location of the Sgr lead- 
ing tail. Note that this part of the stream distributes 
nearly perpendicularly to the Galactic disk in this direc- 
tion of the sky. The panels of Fig. [5] show two different 
3-dimensional projections (top and bottom row) of our 
Sgr model computed for oblate qh = 0.8 (left column) 
and prolate qh — 1.4 halos (right column). They visual- 
ize the simulated Sgr dwarf and its tidally stripped tails 
in rainbow color map, with the color encoding the col- 
umn density of simulation particles (red being the high- 
est, blue the lowest), which are over-plotted in the panels 
as black dots. Note that the core of Sgr dwarf corre- 
sponds to the darkest region owing to the high density of 
N-body particles that are still gravitationally bound to 
the dwarf's remnants. The coordinate system origin is 
located at the Galactic center with the red axis (X) point- 
ing towards the Sun (yellow circle) and the blue axis (Z) 
pointing towards the North Galactic pole. The region in 
which Juric et al. (2006) detect the VOD is marked by 
the white wire-frame box, whereas the location of VSS as 
reported by Duffau et al. (2006; see Sec. 1) is shown by 
the blue wire frame. This coincidence points to a plau- 
sible identification of VOD with the leading arm of the 
Sgr tidal stream or, at the very least, to likely significant 
pollution of VOD as seen in the SDSS density maps by 
the Sgr tidal debris. 

Taken together, the analysis of figures [I] [21 SI and [5] 
clearly suggests that the Virgo over-density may corre- 
spond to the detection of the Sgr stream's leading arm 
falling onto the Milky Way plane. Furthermore, the ob- 
served morphology is best described by models assuming 
an oblate Milky Way halo. With these results at hand, 
it appears fairly unlikely that the VOD is a tidal stream 
unrelated to the Sgr dwarf and yet it has the predicted 
location and morphology of the Sgr leading arm. 

3.2.3. Radial velocity of the Virgo over-density 

At present there is no kinematic information available 
on the stars of Virgo over-density. Radial velocities are 
necessary in order to establish unambiguously whether 
the Virgo over-density is associated to the Sgr stream. 
Theoretical models predict that, independently of the 

11 This figure is a snapshot of a set of movies included as sup- 
plementary material in the electronic version of this paper and also 
available at http://www.astro.princeton.edu/~mjuric/sgr-virgo/ 
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halo shape, the leading arm of the Sgr stream is present 
in this area of the sky and moves on average with nega- 
tive radial velocities as it falls onto the Milky Way disk. 
However, the exact magnitude of the mean radial velocity 
is a vital discriminator for the halo shape, as prolate and 
oblate halos lead to fairly distinct trends of V r along the 
leading arm (Martmez-Delgado et al. 2004a; Helmi 2004, 
Law et al. 2005). We must remark that at present there 
is no information on the kinematics of the Sgr leading 
arm in this area of the sky, since this part of the stream 
was not covered by the spectroscopic survey of M-giants 
carried out by Majewski et al. (2004). 

There have been radial velocity measurements of RR 
Lyrae stars in the Virgo constellation that have shown 
the presence of a spatial over-density (called the 12^.4 
clump, see Sec. 1) located at D ~ 20 kpc. Duffau et 
al. (2006) find a mean radial velocity in the galactic rest 
frame of (V r ) gsr ~ 99.8 km/s (see Table 1) with a fairly 
low velocity dispersion ( er=17.3 km s _1 ), unexpected for 
the smooth population of RR Lyrae in the Milky Way, 
that suggests a possible tidal stream detection. Accord- 
ing to the theoretical models of the Sgr stream, the posi- 
tive value of {y r ) gsr excludes a membership of these stars 
to the leading arm. However, it is unclear whether the 
12 h .4 clump is not associated to the Sgr stream at all, 
since the oblate halo models (qh = 0.8 in this work) also 
predict the presence of the trailing arm stars in that vol- 
ume of the Milky Way (see Fig. [51 where the blue box 
represents the 12 h .4 clump position) with an averaged 
radial velocity of {V r ) gsr = 90 km/s, fairly similar to that 
reported by Duffau et al. (2006) ((V r ) gsr =+83 or +100 
km/s, depending of their sample selection). Therefore, if 
the 12 h .4 clump proves to belong to the Sgr trailing arm 
it would favor the models with < 1. 

4. DISCUSSION 

Theoretical models of the Sgr stream are not able to 
reproduce all the observational constraints. All models 
agree that the geometry of the tidal stream (i.e its 3D 
spatial distribution) clearly rules out halo models with a 
prolate shape (Martmez-Delgado et al. 2004a, Johnston 
et al. 2005, Law et al. 2005) and that Milky Way halos 
with qh < 1 are preferred at a 3tr level. Paradoxically, 
the models that best explain the stream structure fail to 
reproduce the radial velocity trend along the leading arm 
unless prolate halo models (qh > 1) are invoked (Helmi 
2004). Apparently, the present status of the Sgr stream 
modeling is that of a fish biting its tail. 

It is not clear how to reconcile observations and the- 
ory. Possible explanations for this mismatch are: (i) the 
Milky Way models are over-simplistic (although similar 
models successfully reproduce other tidal streams such 
as Monoceros, Penarrubia et al 2005, and the Magellanic 
stream, Ruzicka, Palous & Theis 2006, both streams be- 
ing best reproduced by models with oblate halos) or (ii) 
the presence of possible systematic errors on the data of 
the Sgr tidal stream used in the computation of these 
simulations (as suggested by Chou et al. 2006). 

In this contribution, we have investigated the origin 
of the Virgo over-density discovered from the analysis of 
the SDSS data. The VOD is the largest clump of tidal 
debris ever detected in the outer halo of the Milky Way. 
We find that the position, apparent angular size, stellar 
density and three dimensional geometry match the the- 



oretical predictions from the present Sgr stream models. 
According to these models, the Virgo over-density cor- 
responds to Sgr leading arm falling onto the Milky Way 
disk. In agreement with other portions of the Sgr stream, 
the spatial distribution of the VOD is best reproduced by 
models with oblate halos < 1. However, we must re- 
mark that a unambiguous association between the VOD 
and the Sgr stream cannot be established in absence of 
stellar kinematics. 

If the Virgo over-density proves to belong to the Sgr 
stream it might represent the key piece to solve the para- 
dox described above, as we have for the first time the 
opportunity to measure the radial velocity of millions of 
stream stars just in the region where observations and 
theoretical models do not match. 

Additionally, the detection of a new portion of the 
stream's leading arm provides strong constraints on fu- 
ture theoretical models and will help to address whether 
the Sgr stream passes through the solar neighbourhood 
as it crosses the Milky Way plane (as predicted by Kundu 
et al. 2002, Majewski et al. 2003; Martmez-Delgado et 
al. 2004). In that case, the Sgr stream would repre- 
sent an important target for the on-going RAVE sur- 
vey (Steinmetz et al. 2006), whose goal is to measure 
the radial velocity of around 10 6 stars in both hemi- 
spheres. Thanks to the large sample of radial velocities, if 
the Milky Way crossing occurs at a heliocentric distance 
D < 3 kpc, its presence will be detected in form of an 
excess of stars with v r < in the Northern hemisphere 
(b > 0) and with v r > for b < with respect to the 
background population of Milky Way stars (Penarrubia, 
Navarro & Helmi, in preparation). The presence of a 
tidal stream in the solar neighbourhood also represents 
an excellent target for Weakly Interacting Massive Parti- 
cles (WIMPs) direct detection experiments, since a sub- 
structure of dark matter in the proximity of the Sun 
would yield a "cold" flow of WIMPs thought the detector 
of either of these dark matter experiments (see Savage, 
Freese and Gondolo 2006 and Gondolo et al. 2005 for 
details). This would make the solar vicinity an excellent 
lab for probing the physical properties and nature of the 
dark matter. 
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Fig. 1. — Density maps from SDSS of the North Galactic hemisphere. Position of our positive (solid triangles) and negative 
(open circles) detection of the Virgo stellar over-density given in Table 1, overplotted on a Lambert projection of the density 
map of Virgo given in Juric et al. (2005). This map was generated by selecting MS stars with b < Odeg, 0.2 < g — r < 0.3 
and 20 < r < 21. Concentric circles correspond to constant galactic latitude (the North Galactic pole is in the center; and the 
outermost circle correspond to b=0deg. The high density patch at 1=45, b > 30 is the Milky Way disk protruding into the 
survey region. 



TABLE 1 

Positive and negative detections of Virgo from colour-magnitude diagrams 



R A (,12000) 


Dcc(J2000) 


I 


b 


Detection 


12:24:00.00 


-1:00:00.0 


288.62 


61.11 


Positive 


12:24:00.00 


-0:30:00.0 


288.39 


61.59 


Positive 


10:20:00.00 


40:00:00.0 


180.71 


56.24 


Negative 


13:48:16.00 


52:08:32.0 


103.69 


62.81 


Negative 


15:37:49.00 


69:17:31.0 


104.68 


41.44 


Negative 


12:40:00.00 


18:20:00.0 


285.67 


80.82 


Negative 


13:56:21.00 


-27:10:04.0 


320.28 


33.51 


Negative 


13:07:53.00 


29:21:00.0 


65.42 


85.75 


Negative 


10:49:22.00 


51:03:04.0 


158.57 


56.78 


Negative 


14:07:37.00 


11:50:28.0 


356.10 


66.48 


Negative 


16:11:05.00 


14:57:29.0 


28.76 


42.18 


Negative 


14:44:27.00 


-1:00:00.0 


351.54 


50.88 


Negative 
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Fig. 2. — Galactic X-Z projection of our N-body simulation of the Sgr tidal stream in a Milky Way dark matter halo with 
qn = 0.8. For simplicity, only particles unbound less than 2.5 Gyr ago have been plotted. It displays the Sgr leading arm falling 
onto the Galactic disk, crossing the proximity of the Sun (yellow circle). The particles with a position consistent with that of 
the VOD are marked with red colour. 
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Fig. 3. — Density maps of the Sgr leading tail stars with heliocentric distance d<i5fc pc from our N-body simulation of the Sgr 
tidal stream in a Milky Way dark matter halo with q h = 0.8 (top panel) andq h =1.4 (bottom panel). 
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t = 10 kpc (0.100 < r-i < 0.150) t = 10 kpc (0.10Q < r-i < 0.150) i = 10 kpc (0 100 c r-i < 0.150) 
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Fig. 4. — Top row: comparison of Galactic X-Y projections of the Virgo over-density obtained from SDSS data (left column), 
with two Sgr stream models, qu = 0.8 (center) and qn = 1-4 (right column). In the left column, we plot the residual stellar 
density after subtraction of the Galactic stellar density model in a thin Az ~ 0.5kpc slice at Z = 10 kpc above the Galactic 
plane. The color coding represents the stellar overdensity of 0.1 < r — i < 0.15 stars in stars pc -3 . The colored contours on the 
center and right rows are the Gaussian-smoothed (a — 0.5 kpc) particles from the Sgr simulation (in arbitrary density units). 
To account for the observational errors, the simulation particles were additionally scattered in the radial direction, in accordance 
with the measured SDSS photometric scatter and the intrinsic scatter of the photometric paralax relation (Fig 4., Juric et al. 
2006). Bottom row: comparison of stellar number density in SDSS slice of Virgo over-density with the Sgr simulation. The Z 
axis denotes the distance from the Galactic plane. The R coordinate is the distance from the Galactic center in the <j> — 30° 
direction counter-clockwise along the axis of Galactic rotation. The color coding and the symbols are the same as in the top 
row. 
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Fig. 5. — Two different 3- dimensional projections (top and bottom rows)of the Sagittarius tidal stream models computed for 
an oblate halo qn = 0.8 (left column) and a prolate qn = 1.4 halo (right column). The column density of simulation particles 
is represented by rainbow color-map, with the red being the areas of highest, and blue the lowest density. Over-plotted on the 
column density map as black dots are the simulation particles themselves. The core of Sgr dwarf is clearly visible as a dense, 
dark region, due to the high density of simulation particles which are still gravitationally bound to the dwarf's core. The origin 
of the coordinate system axes is the Galactic center, with the red axis (X) pointing towards the Sun (yellow circle), the blue axis 
(Z) pointing towards the north Galactic pole. The region in which Juric et al. (2006) detect the Virgo overdensity is marked by 
the white wire-frame box. The blue wire-frame marks the location of Virgo Stellar Stream as reported by Duffau et al. (2006). 



